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Abstract We review the effects of macromolecular crowding on the folding of RNA 

S ' 

• ' by considering the simplest scenario when excluded volume interactions between 

O 
f^ • crowding particles and RNA dominate. Using human telomerase enzyme as an 



example, we discuss how crowding can alter the equilibrium between pseudoknot 
and hairpin states of the same RNA molecule — a key aspect of crowder-RNA 
interactions. We summarize data showing that the crowding effect is significant 
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only if the size of the spherical crowding particle is smaller than the radius of 
gyration of the RNA in the absence of crowding particles. The implication for 
function of the wild type and mutants of human telomerase is outlined by using a 
relationship between enzyme activity and its conformational equilibrium. In addi- 
tion, we discuss the interplay between macromolecular crowding and ionic strength 
of the RNA buffer. Finally, we briefly review recent experiments which illustrate 
the connection between excluded volume due to macromolecular crowding and the 
thermodynamics of RNA folding. 

Keywords Crowding ■ Excluded volume • RNA ■ Telomerase ■ Enzyme activity 
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1 Introduction 



The cytosol is crowded and replete with macromolecules such as ribosomes, hpids, 
proteins, and RNA. Estimates show that the volume fraction {4>) of these macro- 
molecules, collectively referred to as crowding agents, can exceed 0.2. Spontaneous 
folding of nascent proteins and RNA in the crowded environment can be differ- 
ent from in vitro experiments that are typically conducted under infinite dilu- 
tion. Minton, who has made pioneering contributions in elucidating the impor- 
tance of crowding in biophysics, was the first to recognize that the thermodynam- 
ics of f olding, association, and biochemical reactions can be altered by cr owding 



agents (JMinton , 



1981 



2001 



2005 



Hall and Minton, 



2003: 



Zhou et al 



2008 ). More 



recently there has been much interest in the study of crowding effects on folding 



and function of proteins (| Cheung et al , 



Elcock . 



2005 



Homouz et al . 



2008; 



Dhar et al . 



2010 : 



2010|) , which was inspired by the insightful studies initiated by Minton. 



Describing even the simplest process, namely, the transition between folded 
and unfolded states of proteins and RNA under crowded conditions is complicated 
because the nature of the effective interactions between the crowding agents and 
the polypeptide or polynucleotide chains is not fully understood. However, to a first 
approximation, the dominant effect of crowding agents is to exclude the molecule 
of interest from the volume occupied by crowders. If excluded volume interactions 
dominate (an assumption that has to be tested before the theory can be applied to 
analyze experiments) , then the stability of the folded state of the protein or RNA 
is enhanced, compared to the case when (^ = 0. In this case, the loss in entropy 
of the folded state due to crowding is much less than that of the unfolded state, 
resulting in the stabilization of the folded state. The entropic stabilization of the 
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folded state (jCheung et al , 



2005 : 



Minton 



2005|) in the presence of crowding agents 



has been firmly established theoretically in a number of studies. 

Most of the studies on the effects of crowding on self-assembly process have 
been on the folding of proteins. Recently, it has been recognized that crowd' 



ing agents could have a significant impact on the folding of RNA (jPincus et al , 



2008 : 



Kilburn et al . 



201C; 



Denesvuk and Thirumalai , 



20111) . Simple theoretical ar- 



guments and coarse-grained simulations wer e used to show th at crowding can 



2008 ). However, RNA 



modestly stabilize RNA secondary structures (jPincus et all 
requires counterions (Mg^"*", for example) for tertiary folding. Thus, the effect of 
macromolecular crowding on tertiary structures of RNA may be complicated. Us- 
ing small angle X-ray scattering measurements it has been shown that, in the 
presence of polyethylene glycol ( PEG), the 195 nucl eotide Azoarcus ribozyme is 



more compact relative to (f> = (JKilburn et al 



20101) . It was concluded that ex- 



cluded volume effects play a dominant role in the compaction of RNA in low 
molecular weight PEG. Interestingly, the transition to the folded state o ccurs at 



20101) . Even if 



a lower Mg^"*" concentration in the presence of PEG ( Kilburn et al 

excluded volume interactions largely determine the stability of the folded states of 

RNA, a number of variables besides 4>, such as size and shape of crowding agents, 

also contribute to the stability of RNA in the presence of inert crowding agents. 

Thus, a systematic study of the infiuence of macromolecular crowding on RNA is 

required. 



It is known that, in contrast to proteins (JGuo et al 



1992), the s tability gap sep- 



arati ng the native state and low energy excitations in RNA is small (JThirumalai and Hveon , 



2005r) . which implies that external factors (crowding or force) can modulate the 



stabilities and functions of RNA. For example, riboswitches undergo a transition 
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betwe en two distinct conformation s that have a profound influence on their func- 



tions (JMontange and Batev . 



20081) . In this review article, we outline essential as- 
pects of crowder-RNA interactions, with the particular emphasis on how crowding 
in the cellular environment may alter conformational equilibria related to enzy- 
matic function. As a biologically relevant example, we discuss crowding effects on 
the transition between the hairpin (HP) and pseudoknot (PK) conf ormations (Fig 



ure[T] ) in the pseudoknot domain of human telomerase RNA, hTR (JTheimer et al . 



20031) . The pseudoknot domain is conserved i n different organisms and its ac- 



tivity is closely linked to chromosome stability (JBlasco , 



2003 : 



Chen and Greider . 



20041) . However, the precise role of the PK and HP conformations of the pseudo- 
knot domain in the context of telomerase activity is not known. Mutations that 
either increase or decrease t he stability of the PK conformation resuh in a reduc- 



tion in telomerase activity (jComoUi et a. 



2002 : 



Theimer et a' 



20051) . Therefore, 



it is important to compare the impact of physical factors, such as macromolecular 
crowding, with that of naturally occurring chemical mutations. Our discussion of 
crowding presented in this review is original in two ways. First, we consider the 
effects of crowding on the conformational equilibrium between two folded states. 
Traditionally, discussions of crowding effects on biomolecules, mostly proteins, 
concern with the stability of the folded (active) state with respect to the unfolded 
(inactive) state. In addition, we establish a quantitative and novel connection be- 
tween the magnitude of the crowding effect and activity of wild-type and mutant 
enzymes. 



Natalia A. Denesyuk, D. Thirumalai 



2 Computational models of RNA and crowders 

Entropic stabilization, the mechanism by which macromolecular crowding stabi- 
lizes folded states of proteins and RNA in the excluded volume limit, is similar to 
spatial confinement, in a sense that crowders confine the biomolecule of interest 
to the interstitial space. Although analytical results exist for polymers confined to 
cavities with simple geometries, it is difficult to obtain accurate estimates for non- 
trivial confinement geometry associated with macromolecular crowding. Further- 
more, unfolded RNA is a highly nonideal polymer whose conformational ensemble 
will be determined, among other things, by stacking propensities between adjacent 
bases and by the ionic strength of the buffer. Estimating how the conformational 
space of such a polymer will be reduced even in a simple geometrical confinement 
is a challenging task. 

On the other hand, simulations have proven to be a useful tool in assess- 
ing the magnitude of entropic stabilization in response to varying external condi- 
tions. Coarse-grained models have been particularly effective in the studies of 
RNA folding, since they do n ot suffer from computational complexity associ' 



ated with all-atom force fields ( 



Cho et al . 



2005 : 



Feng et al 



Lin and Thirumalai. 



20111 



2008; 



Whitford et a" 



Denesvuk and Thirumalai , 



2009 : 



20111) . As a specific 



example of this class of techniques, we will be discussing results obtained with a 
widely appreciated coarse-grained representation of RNA, where each nucleotide 
i s modeled by three interactions sites (TIS) — a phosphate, a sugar and a bas e 



( Hveon and Thirumalai , 



2005 : 



Cho et al . 



2008 



Denesvuk and Thirumalai , 



2011) 



We developed a force field in conjunction with the TIS model, which originally 
included stacking and hydrogen bond interactions as essential components in the 
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stability of RNA st ructures and was thermodynami cally accurate in the limit of 



high ionic strength (JDenesvuk and Thirumalai , 



20111 ) . The quantitative agreement 



of thermodynamic predictions of this force field with experiments is illustrated in 
Figure [3^. Subsequently, we added an electrostatic component to the TIS model 
to describe the RNA thermodynamics at different ionic concentrations (Denesyuk 
and Thirumalai, to appear elsewhere). The description of electrostatic effects in 
the model is based on Manning's concept of counterion condensation, which posits 
that counterions condense onto the RNA molecule and reduce the charge of phos- 
phate groups from — e to — Qe, where Q < 1 and e is the proton charge. The 
uncondensed mobile ions are treated using the Debye-Hiickel theory. Therefore, 
the electrostatic energy of RNA in simulation is computed using 



^Vv- 



exp(-|r; -rjl/A) 



^^ = ^^1. "\. ' /'^^ (1) 



2s ^^ 

1,3 



where |rj — r ^ | is the distance between two phosphates i and j, e is the dielectric 
constant of water and A is the Debye-Hiickel screening length. We showed that, if 
the reduced phosphate charge Q was taken to be 



Q = r, (2) 



where l^ is the Bjerrum length and b = 4.4 A, the measured thermodynamics 
of a variety of RNA sequences could be reproduced well over a wide range of 
temperatures and monovalent salt concentrations (Denesyuk and Thirumalai, to 
be published). 

A generalized Lennard-Jones potential, introduced by 



Denesvuk and Thirumalai 



(J201lf ). has been successfully employed to model interactions of RNA with the 



Natalia A. Denesyuk, D. Thirumalai 



spherical crowders of arbitrary size, 

Do V^ o/ Do ^6 



^0 



( ^0 Y^ 2( ^° V I 1 

\r + Do-Dj \r + Do-Dj 



r <D, 



Ui^jir) = 0, r > D, (3) 

where r is the distance between the centers of mass of two interacting particles, 
Do is the effective penetration depth of the interaction, Ri is the radius of an 
RNA coarse-grained bead, re is the radius of a crowder, and D = Ri + re- The 
ratio 2Ri/Do in Eq. ^ is used to scale the interaction strength e in proportion to 
the surface contact area. This potential correctly accounts for nonspecific surface 
interactions between spherical crowders representing large macromolecules and 
individual segments of the coarse-grained RNA. 

Assessing the magnitude of crowding effect from simulations requires an ac- 
curate technique for calculating the folding free energies. To this end, a thermo- 
dynamic method has been proposed (Denesyuk and Thirumalai, to appear else- 
where), which does not require a structural order parameter to define folded and 
unfolded ensembles. To summarize, we perform a series of simulations at different 
temperatures in the range from Ti to T2, where Ti and T2 are on the order of 
the lowest and highest temperatures used in thermodynamic simulations or mea- 
surements. Using statistical mechanics techniques, we compute from these simu- 
lations the free energy of the molecule, G, as a function of temperature, T (Figure 
[2]) . If the population of the unfolded state is negligible at Ti , the free energy of 
the folded state at Ti, Gf(Ti), equals the computed free energy G(Ti). Similarly, 
Gu{T2) = G(T2), if the population of the folded ensemble is statistically insignifi- 
cant at the highest temperature T2. Assuming similar equalities for the enthalpies 
and heat capacities of the folded and unfolded state, we can use thermodynamic 
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relationships between the free energies, enthalpies and heat capacities to extrap- 
olate Gf (Ti) and Gu{T2) to intermediate temperatures, at which both folded and 
unfolded states are populated (asymptotic lines in Figure [2]). For any T between 
Ti and T2, geometric definition of the stability of the folded state with respect to 
the unfolded state, AG, is given in Figure [2] 



3 Crowding effect is negligible for large crowders 



The high density of macromolecules in the cell (volume fractions 4> Ki 0.2 — 0.4) re- 
duces the space available for conformational fluctuations. Therefore, macromolec- 
ular crowding should result in a shift in the thermodynamic equilibrium between 
the HP and PK states of the hTR pseudoknot domain towards the more com- 
pact PK. To assess the extent to which PK i s favor ed at </) 7^ 0, we discuss the 



simulation results (JDenesvuk and Thirumalai , 



20111 ) for the HP and PK states 



of the modified pseudoknot domain, A\]177 (Figure [T|). Th e molecular construc t 



2005 1 



ziU177 has been examined experimentally in vitro at = (jTheimer et al 
and the atomistic structures of its HP and PK conformations are available from 
the Protein Data Bank, codes 1NA2 and 2K96, respectively. Under native con- 
ditions, the RNA sequence in Figure [1] will predominantly populate the PK con- 
formation. To obtain adequate statistics on both conformations, two independent 
sets of sim ulations were carried out, each m odeling a limited subset of hydro- 



gen bonds (JDenesvuk and Thirumalai , 



201lh . In the first set of simulations, the 



RNA sequence could form only those hydrogen bonds that are found in the NMR 
structure of the PK. Similarly, in another set of simulations, only hydrogen bonds 
from the NMR structure of the HP were included. In the HP simulations strand 
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C166-A184 remained unbound in the folded state, because no hydrogen bonds 
could form between this strand and the remainder of the molecule. In this way, 
the interconversion between the PK and HP was eliminated and the HP struc- 
ture could be sampled exhaustively. The stabilities of the PK and HP structures, 
AGpK = GpK,f — GpK.u and AGup = GHP,f — Ghp.u, were obtained from the 
simulations using the technique illustrated in Figure [2] Since the unfolded state is 
effectively defined as a high-temperature state (in which all hydrogen bonds are 
broken) and the RNA sequence is the same in the PK and HP simulations, it was 
assumed that Gpk.u = Ghp,u- Therefore, the stability of the PK with respect to 
HP, GpK.f — Gjjp.f , was computed as Z\Gpk — ^^Ghp, without the need for explicit 
simulations of the PK-HP interconversion. 

In order to illustrate the essential aspects of crowding effects on RNA, we 
consider spherical particles with radius vq- For monodisperse particles, the vol- 
ume fraction is (f> = 47rr^p/3, where p is the number density. Thus, tj) can be 
changed by increasing or decreasing p or by altering the size of the crowding 
particles. In this review, we fix = 0.3 and examin e the consequences of c hang- 



ing re- Based on general th eoretical considerations (| Asakura and Oosawa , 



Shaw and Thirumalai, 



1958 : 



1991 ) it can be shown that, in the colloid limit re > Rq, 
the crowding agents would have negligible effects on RNA stability. Here, _Rq is 
the size of RNA in the absence of the crowding agent. It is only in the opposite 
polymer limit, re < Rq, that the crowding particles would affect RNA stability. 
We therefore expect that the magnitude of the crowding effect should depend on 
the ratio tc/Rq- 

In Figure[3^ we show the HP melting profile, taken to be the negative derivate 
of the number of intact base pairs A'bp with respect to k^T, for the crowder radius 
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re = 26 A ( Denesvuk and Thirumalai , 



201lf) . Such crowders are larger than the 



radius of gyration of strand G93-C121 in the unfolded state, Rq = 20 A (Figure 
[TJ;) . As discussed above, large crowders have minimal effects on the melting of the 
HP even at = 0.3. For a fixed 4>, the average distance between two spherical 
crowders will increase with the crowder size. If the unfolded hairpin can easily 
fit in the interstitial space, the folding/unfolding transition will not be affected 
significantly by the presence of crowders. For (j> = 0.3 and the crowder radius 
re = 26 A, which is only slightly larger than i?Q, the increase AT in the melting 
temperature is 1.5 °C for stem 1 of the HP and is negligible for stem 2 (Figure 
|3^). Further increase in re results in AT f» for both stems (data not shown). 

Figure[3^ also shows the melting profile of the HP in a ternary mixture of crow- 
ders, containing volume fractions </> = 0.11, 0.11 and 0.08 of particle s with re = 104 



A, 52 A and 26 A, respectively ( Denesvuk and Thirumalai , 



201l|). The sizes and 



volume fractions of individual components in the model mixture correspond to 
the ribosome, large enzymatic complexes and relatively small individual proteins, 
found in E. coli. Because all the values of re in the E. coli mixture are larger than 
Rq , we expect only small changes in the melting profile of the HP (Figure [3^) . 
For the total volume fraction of 0.3, the melting temperature of the HP stem 1 
increases only by 2 °C with respect to 4> = (Figure |3^). Interestingly, the effect 
of the E. coll mixture is similar in magnitude to that of a monodisperse suspension 
with re = 26 A and tj) = 0.3. In contrast, a monodisperse suspension with re = 26 
A and (j) = 0.08, which is equivalent to the smallest particle component in the 
mixture, has negligible effect on the melting of the HP (Figure [3^). 

To summarize, the crowding effect of polydisperse mixtures is largely the effect 
of the smallest particle component, but taken at the total volume fraction of the 
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mixture. As we discuss in the next section, tlie excess stability of the folded state 
due to crowding decreases nonlinear ly with the radius of the crowding particle re- 
We therefore propose that, for crowding in the cellular environment, the main role 
of large macromolecules will be to increase the effective volume fraction of the 
relatively small macromolecules. 



4 Role of crowder size in the PK-HP equilibrium 

The sensitivity of the crowding effect to the relative sizes of RNA and crowders 
is at the basis of the equilibrium shift in the hTR pseudoknot domain. Figure [3}3 
shows the change in stability of the HP and PK at 37 °C induced by monodisperse 
crowders for different crowder radii re (0 = 0.3). As anticipated by arguments 
given above, the magnitude of the excess stability Z\G(0.3) — AG{0) is small if 
tc/Rg > 1 and increases sharply for tc/Rg < 1- Note that the crowding effect is 
larger for the PK for all values of re (Figure |3)d) , indicating an equilibrium shift 
towards this conformation. The discussed change in the PK-HP relative stability 
is of entropic origin. The unfolded ensembles of the PK and HP are equivalent, 
as discussed above, and will therefore be depleted by macromolecular crowding to 
a similar degree. The compact folded PK is, to a first approximation, unaffected 
by crowding. However, the folded HP contains an unbound strand Cf66-Af84 
(Figure[T|), whose loose conformations will be restricted in a crowded environment. 
Therefore, the excess stability of the PK with respect to HP is due to a partial 
suppression of the HP folded ensemble by macromolecular crowding. 

The crowder radius re = 12 A corresponds to the size of an average protein 
m VIVO. For (j) = 0.3 and re = 12 A, we have ziG'pK(0.3) - Z\Gpk(0) = -2.4 
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kcal/mol and Z1Ghp(0.3) — Z\G'hp(0) = —1.0 kcal/mol, which amounts to the 
relative stabihzation of tlie PK conformation by —1.4 kcal/mol (Figure[2b)- Below 
we analyze this value in the context of standard changes in the PK stability caused 
by mutations. 



5 Implications for function 



As mentioned in the Introduction, changes in the relative stability of the HP and 
PK conformations of the hTR pseudoknot domain compromise the enzyme activ- 
ity. The estimate of the crowding effect in a typical cellular environment, A AG = 
-1.4 kcal/mol, allows us to assess the extent to which macro molecules could regu- 



late telomeras e activ ity. Experimental data on hTR mutants (|Comolli et al . 



Theimer et al . 



2002: 



2005 ) clearly demonstrate that enzyme activity decreases when 
plotted as a function of \AAG*\ = \AGp^{0) - ZiGpK(0)|, where AGpj^{0) and 
Z\Gpk(0) are the stabilities of mutant and wild-type pseudoknots at = (Figure 
H]). The majority of mutations destabilize the PK, A AG* > (black squares in Fig- 
ureSI and only two mutants have AAG* < (red stars in FigureBl). For destabiliz' 



i ng m utants the reduction in activity, a, was shown by 



Denesvuk and Thirumalai 



(J2011I ) to follow the exponential dependence, a = exp(— 0.6Z\Z\G*) (thick curve 
in Figure 131). The naturally occurring destabilizing mutations DKC and C116U 
have b een linked to d i seases dyskeratosis congen ita and aplastic anemia, respec- 



tively (jVulliamv et al . 



2002 



Fogartv et al . 



20031 ). The DKC and the stabihzing 



ZiU177 mutations have been studied in vivo (green symbols in Figure |4|), as well 
as in vitro. In both cases, mutant telomerase m vivo was found to be significantly 
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less active than the corresponding construct in vitro, suggesting that a number of 
factors deterinine the activity of telomerase in vivo. 

Although macromolecular crowding enhances the stability of the PK state, the 
crowding effect {A AG = — 1.4 kcal/mol) is less than the stability changes caused 
by mutations. In Figure |4] the grey area marks the domain of potential mutants 
with A AG* > 0, whose activity may be completely restored by macromolecular 
crowding. All experimentally studied mutants fall outside the marked domain, in- 
cluding the two disease related mutants DKC and C116U. Nevertheless, due to the 
strong dependence of enzyme activity on A AG* , the effect of crowding on telom- 
erase function may be significant. We estimate that the activity of telomerase can 
be up- or down-regulated by more than two-fold in response to density fluctua- 
tions in its immediate environment. Furthermore, due to the expected dynamical 
heterogeneities in cells, there will be variations in enzyine activity in different cell 
regions. 

6 Crowding effects on RNA at different ionic strengths 



Entropic stabilization mechanism ( Cheung et al . 



20051) implies that crowding in- 



creases the stability of the folded state by reducing the population of expanded 
conformations in the unfolded state. Therefore, we expect that the magnitude of 
the crowding effect will be sensitive to the ionic strength of the RNA buffer, since 
the latter determines the size of the unfolded RNA. The quantitative discussion 
above assumed the limit of high ionic strength. As the buffer ionic concentration 
c is lowered, the screening of the negative charge on the RNA sugar-phosphate 
backbone becomes less efhcient, which in turn increases the mean radius of gyra- 
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tion of conformations in the unfolded state. The function -Rg(c) for the unfolded 
PK is shown in Figure [S^ in the absence (black diamonds) and presence (green 
circles) of crowding. The same general trend is observed in both cases, with the Rq 
values being consistently smaller when crowders are present for the entire range of 
c. In accord with our predictions, the crowder-induced stabilization of the folded 
PK becomes more significant at low ionic strengths (red squares in Figure [5^) . 
Interestingly, the stabilization effect increases rapidly upon lowering c from 1 M 
to 0.1 M, but shows little change when c is lowered further to below 0.1 M. The 
underlying reason for such behavior can be traced to the probability distributions 
p{Rq) in the unfolded state (Figure[5}3) . For a given crowder solution, we can iden- 
tify a typical size of the cavity which will be free of any crowders. If the radius of 
gyration of RNA conformations is such that they fit into the cavity, these confor- 
mations will not be perturbed by crowding. On the other hand, the population of 
conformations with Rq larger than the typical cavity size will be significantly de- 
pleted by crowding. For (f> = 0.3 and re = 12 A, we can infer the size of a standard 
empty cavity from the distributions p{Rq) at high ionic strength (c = 1 M, Figure 
[5]d). Note that p(-Rg) decreases for Rq > 20 A when crowders are present (green 
solid line in Figure[5]3), but p{Rq) increases with Rq around 20 A in the absence 
of crowders (black solid line in Figure [SJd). This indicates that the crowders signif- 
icantly perturb the unfolded conformations with _Rg larger than 20 A, which can 
serve as an upper estimate of the smallest RNA size affected by crowders. When 
c decreases, the distribution p{Rq) shifts to larger _Rgi increasing the fraction of 
the unfolded conformations affected by crowders. At c = 0.1 M, all statistically 
significant values of Rq in the unfolded state fall within the range Rq > 20 A, 
both for cj> = (black diamonds in Figure \Ejp) and for (p = 0.3 (green circles in 
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Figure[5}3), so that the entire distribution p{Rg) is depleted due to crowding. This 
explains why the crowding-induced stabilization is almost constant below 0.1 M, 
even if the average Rq continues to increase rapidly all the way to 0.01 M (symbols 
in Figure [5^). 



7 RNA becomes compact as </> increases 



The reduction in conformational space accessible to RNA should increase with 4>, 
for a fixed re. Thus, we expect that Rq s hould decrease as (f > increases. This is 



2010l) . which show that 



precisely what is observed in experiments (jKilburn et all 
at all concentrations of Mg ^ the Azoarcus ribozyme becomes more compact as the 
volume fraction of the crowding agent (PEG) increases (Figure [6l). Interestingly, 
the midpoint of the folding transition Cm — the concentration of Mg ^ at which 
the folded and unfolded states of the Azoarcus ribozyme have equal populations — 
also decreases as 4> increases (Figure [6]). This fi nding can be readil y explained in 



20051 ) and suggests 



terms of the entropic stabilization mechanism ( Cheung et all 
that, to a first approximation, PEG behaves as an inert hard sphere crowding 
agent. Based on our considerations from the previous section, we predict that the 
shift in Cm due to crowding will also depend on the concentration of monovalent 
counterions in the RNA buffer, a prediction that is amenable to experimental test. 
In addition, it would be of interest to perform experiments at a fixed but varying 
re, which can be changed by decreasing or increasing the molecular weight of PEG. 
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8 Conclusions 

The phenomena discussed here illustrate just one aspect of crowding effects on 
RNA. There are other physical and chemical factors that could determine how 
crowding modulates RNA stability, and hence its function. The interplay between 
electrostatic interactions, hydration of phosphate groups and crowding particles, 
as well as the shape of crowding particles, could potentially present a much more 
nuanced picture than the simplest case considered here. Even though these effects 
are important and warrant further study, the analysis presented here is important 
in establishing the maximuin increase in stability that can be realized when ex- 
cluded voluine interactions dominate. From this perspective, our conclusion that 
the crowding effects cannot fully restore telomerase activity in vivo is robust. In 
future work, it will be important to connect the consequences of crowding effects 
on RNA folding with the RNA activity under cellular conditions. 

Acknowledgements Wc are grateful to Sarah Woodson for providing Figure 6. Our work 
was supported by a grant from the National Science Foundation (CHE 09-10433). 

Conflict of Interest: None 



References 

Asakura S, Oosawa F (1958) Interaction between particles suspended in solutions 

of macromolecules. J Polym Sci 33(126):183-192 
Blasco M (2003) Telomeres and cancer: a tale with many endings. Curr Opin Genet 

Dev USA 13(l):70-76 



18 Natalia A. Denesyuk, D. Thirumalai 

Chen J, Greider C (2004) An emerging consensus for telomerase RNA structure. 

Proc Natl Acad Sci USA 101(41):14,683-14,684 
Cheung M, Klimov D, Thirumalai D (2005) Molecular crowding enhances native 

state stability and refolding rates of globular proteins. Proc Natl Acad Sci USA 

102(13):4753-4758 
Cho SS, Pincus DL, Thirumalai D (2009) Assembly mechanisms of RNA pseu- 

doknots are determined by the stabilities of constituent secondary structures. 

Proc Natl Acad Sci USA 106(41): 17,349-17,354 
ComoUi L, Smirnov I, Xu L, Blackburn E, James T (2002) A molecular switch 

underlies a human telomerase disease. Proc Natl Acad Sci USA 99(26): 16, 998- 

17,003 
Denesyuk NA, Thirumalai D (2011) Crowding promotes the switch from hair- 
pin to pseudoknot conformation in human telomerase RNA. J Am Chem Soc 

133(31):11,858-11,861 
Dhar A, Samiotakis A, Ebbinghaus S, Nienhaus L, Homouz D, Gruebele M, 

Cheung MS (2010) Structure, function, and folding of phosphoglycerate kinase 

are strongly perturbed by macromolecular crowding. Proc Natl Acad Sci USA 

107(41):17,586-17,591 
Elcock AH (2010) Models of macromolecular crowding effects and the need for 

quantitative comparisons with experiment. Curr Opin Struct Biol 20(2):196— 

206 
Feng J, Walter NG, Brooks CL III (2011) Cooperative and directional folding of 

the preq(l) riboswitch aptamer domain. J Am Chem Soc 133(12):4196-4199 
Fogarty P, Yamaguchi H, Wiestner A, Baerlocher G, Sloand E, Zeng W, Read 

E, Lansdorp P, Young N (2003) Late presentation of dyskeratosis congenita 



Entropic stabilization of the folded states of RNA due to macromolecular crowding 19 

as apparently acquired aplastic anaemia due to mutations in telomerase RNA. 

Lancet 362(9396):1628-1630 
Guo Z, Thirumalai D, Honeycutt J (1992) Folding kinetics of proteins — a model 

study. J Chem Phys 97(l):525-535 
Hall D, Minton AP (2003) Macromolecular crowding: qualitative and semiquan- 
titative successes, quantitative challenges. BBA-Proteins Proteom 1649(2); 127- 

139 
Homouz D, Perham M, Samiotakis A, Cheung MS, Wittung-Stafshede P (2008) 

Crowded, cell-like environment induces shape changes in aspherical protein. Proc 

Natl Acad Sci USA 105(33):11,754-11,759 
Hyeon C, Thirumalai D (2005) Mechanical unfolding of RNA hairpins. Proc Natl 

Acad Sci USA 102(19):6789-6794 
Kilburn D, Roh JH, Guo L, Briber R, Woodson SA (2010) Molecular crowding 

stabilizes folded RNA structure by the excluded volume effect. J Am Chem Soc 

132(25) :8690-8696 
Lin JC, Thirumalai D (2008) Relative stability of helices determines the folding 

landscape of adenine riboswitch aptamers. J Am Chem Soc 130(43):14,080- 

14,081 
Minton A (2001) The influence of macromolecular crowding and macromolecu- 
lar confinement on biochemical reactions in physiological media. J Biol Chem 

276(14):10,577-10,580 
Minton AP (1981) Excluded volume as a determinant of macromolecular structure 

and reactivity. Biopolymers 20(10) :2093-2120 
Minton AP (2005) Models for excluded volume interaction between an unfolded 

protein and rigid macromolecular cosolutes: Macromolecular crowding and pro- 



20 Natalia A. Denesyuk, D. Thirumalai 

tein stability revisited. Biophys J 88(2):971-985 
Montange RK, Batey RT (2008) Riboswitches: Emerging themes in RNA structure 

and function. Ann Rev Biophys 37:117-133 
Pincus D, Hyeon C, Thirumalai D (2008) Effects of trimethylamine n-oxide 

(tmao) and crowding agents on the stability of RNA hairpins. J Am Chem 

Soc 130(23) :7364-7372 
Shaw M, Thirumalai D (1991) Free polymer in a colloidal solution. Phys Rev A 

44(8):R4797-R4800 
Theimer C, Finger L, Trantirek L, Feigon J (2003) Mutations linked to dyskeratosis 

congenita cause changes in the structural equilibrium in telomerase RNA. Proc 

Natl Acad Sci USA 100(2) :449-454 
Theimer C, Blois C, Feigon J (2005) Structure of the human telomerase RNA 

pseudoknot reveals conserved tertiary interactions essential for function. Mol 

Cell 17(5):671-682 
Thirumalai D, Hyeon C (2005) RNA and protein folding: Common themes and 

variations. Biochemistry 44(13):4957-4970 
VuUiamy T, Marrone A, Dokal I, Mason P (2002) Association between aplastic 

anaemia and mutations in telomerase RNA. Lancet 359(9324) :2168-2170 
Whitford PC, Schug A, Saunders J, Hennelly SP, Onuchic JN, Sanbonmatsu KY 

(2009) Nonlocal helix formation is key to understanding s-adenosylmethionine-1 

riboswitch function. Biophys J 96(2):L7-L9 
Zhou HX, Rivas G, Minton AP (2008) Macromolecular crowding and confinement: 

Biochemical, biophysical, and potential physiological consequences. Ann Rev 

Biophys 37:375-397 



Entropic stabilization of the folded states of RNA due to macromolecular crowding 



21 



(a) , ACA, (C) 

A^ ^i/4 Stem 2 ip ,, 



174 



183 
I 



^ACAAAAAAGUCAGCA ^ 

C-166 AAAGUCAGCA" A c.^^ y 3' 

.2l-CCCGACUUUCAGUCG-.07_^^C^l«6,^^^^'-^^Ur.,„, 

3 ,GGGCUGyiJUUUC ^ C 
99 ^ ^ ^ '104 93 'S'^em 1 99 104 



,GGGCUG C 

92 Stem 1 YUuttUQ^ 

QQ ^ U ^ ^1 




Fig. 1 Graphic adopted from 



Denesvuk and Thirumalail 1 2011 ). Secondary and tertiary struc- 



tures of the pseudoknot (PK) and hairpin (HP) conformations for Z\U177. The sequence of 
ziU177 was chosen to match e xperiments, in which nucleotides 122-165 and 177 were deleted 



to facilitate structural studies IITheimer et a] 



2OO5I I . (a) PK secondary structure, (b) PK ter- 



tiary structure, (c) HP secondary structure, (d) HP tertiary structure. NMR structure of the 
HP includes residues G93 to C166 only (PDB code 1NA2). To quantify the effect of crowdcrs 
on the PK-HP equilibrium, we added an unstructured tail A167-A184 to the NMR structure 
(panel c). This inclusion ensures that identical RNA sequences are used in simulations of the 
PK and HP conformations. 
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Fig. 2 Geometrical definition of the stability AG of the folded state. The solid curve illustrates 
the total free energy of a system as a function of temperature, computed in simulation. The free 
energies of individual folded and unfolded ensembles (dashed curves) are obtained as described 
in the text. 
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Fig. 3 (a) Melting profiles of the HP in various crowding environments at 1 M or higher 
monovalent salt concentration. Black solid curve: without crowdcrs. Red dashed curve (the 
model E. coli mixture): = 0.11 of crowdcrs with re = 104 A, if> = 0.11 of crowdcrs with 
re = 52 A and (p = 0.08 of crowdcrs with re = 26 A. Green dashed-dotted curve: 4> = 0.08 of 
crowdcrs with re = 26 A. Blue dotted curve: (p = 0.3 of crowdcrs with re = 26 A. The thick 



Theimer et a] 



orang e curve is the experimental UV data at 200 mM KCl from Figure 2b in 

i 20031 ). divided by 5.53 X 10~^. The two peaks indicate melting of stems 1 and 2 of the HP 



(Figure ^). The peak positions (melting temperatures) and the overall width of the melting 
curve (melting range) servo as a measure of agreement between theory and spectroscopic data, 
(b) Changes in stability (kcal/mol) of the HP and PK at 37 °C due to crowdcrs at </> = 0.3, 
as a function of the crowder radius re • 



24 



Natalia A. Denesyuk, D. Thirumalai 



100-1 



AU177 

l>'=100exp(-3/5x) 



DKCi iC116U 




4 6 8 

|AAG*|, kcal/mol 



Fig. 4 Graphic adopted from 



Denesyuk and Thirumalail ( 20111 ') . Activity of mutant tclomcrasc 



normalized to wild-type activity (100%), as a function of the magnitude of the stability differ- 
ence between mutant and wild-type pseudoknots. Three mutants, Z\U177, DKC and C116U, 
are explicitly marked. The in vitro d ata for destabil i zing ( black squares) and stabilizing (red 



stars) mutations are from Table 2 in 



data for Z\U177 and DKC from 



Comolli et all i 



Theimeretajl 1 20051 '). Green symbols show the in vivo 



2002f ). The experimental data for destabilizing 



mutations is fit to the exponential function (solid curve) . Grey area marks the range of stability 
differences that can be accommodated by crowding. 
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Fig. 5 (a) The radius of gyration of the unfolded PK (right axis), as a function of the mono- 
valent ion concentration c, for </> = (black diamonds) and </) = 0.3, re = 12 A (green circles). 
Red squares show the excess stability of the PK due to crowding, Z\Gpk(0.3) — Z\Gpk(0), 
for re = 12 A (left axis), (b) ProbabiUty distributions p(-Rg) of the radius of gyration of 
the unfolded PK for (^ = and c = 1 M (black soUd line), for (j> = and c = 0.1 M (black 
diamonds), for (j> = 0.3, re = 12 A and c = 1 M (green soUd Hne), and for </) = 0.3, re = 12 
A and c = 0.1 M (green circles). The vertical dashed Hne indicates the smallest size of RNA 
conformations that will be perturbed by crowders with (j> = 0.3 and re = 12 A. 
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Fig. 6 Graphic adopted from 



Kilburn et all 1 2010 ) . Small angle X-ray scattering measurements 



of the radius of gyration of the Azoarcus ribozymc for different concentrations of Mg^+ ions 
and PEG. 



